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Abstract 

Stability of nuclei beyond the drip lines in the presence of an envelop- 
ing gas of nucleons and electrons, as prevailing in the inner crust of a 
neutron star, is studied in the temperature-dependent Thomas- Fermi 
framework. A limiting asymmetry in the isospin space beyond which 
nuclei cannot exist emerges from the calculations. The ambient con- 
ditions like temperature, baryon density and neutrino concentration 
under which these exotic nuclear systems can be formed are studied in 
some detail. 
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I. INTRODUCTION 



In the outer region of the core of a neutron star, at densities ~ 1.5 x 10 14 g/cm 3 
or more, nucleons are distributed uniformly forming a homogeneous system [[[]-§}]. 
At lower densities, in the inner crust of the star, inhomogeneities appear and then 
the positive charges get concentrated in individual nuclei of charge Z which are 
embedded in a sea of neutrons, electrons and possibly neutrinos, the whole system 
being charge neutral. The nuclei arrange themselves in a periodic lattice || due to 
the electrostatic interaction. In the typical conditions prevailing in the neutron star 
interior, the nuclei are in complete thermodynamic equilibrium with the environment 
and are assumed to be in /3-equilibrium. They may well be beyond the neutron drip 
line (THJ known for laboratory nuclei with (N — Z)/A < 0.35. The excess neutrons 
of the very neutron rich nuclei present in the inner neutron star crust, which would 
otherwise decay under terrestrial conditions, are held together in stable equilibrium 
by the pressure exerted by the surrounding neutron sea. 

The properties of the isospin-rich nuclei may be quite different from those of the 
terrestrial nuclei. For example, with increasing density, the nuclei may pass through 
different exotic shapes, namely from spheres to cylinders, slabs, cylindrical holes and 
spherical bubbles . The presence of nonspherical nuclei could affect significantly 
the pinning of vortices and neutrino emission from neutron stars |J. The existence 
of these exotic-shaped nuclei is, however, model-dependent PJIUfl. The external gas 
surrounding the nuclei may also influence the density distribution of these nuclei 
and thus may modify their size [p]]. 

The nuclear equation of state (EOS) plays the pivotal role in determining the 
macroscopic properties of a neutron star, such as its mass, radius and moment 
of inertia [11]]. In this context, the presence of different nuclear species in the 
outer and inner crustal regions provides the starting point for the consideration of 
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some important aspects, such as its superfruid and elastic properties. On the other 
hand, for nuclei immersed in a nucleonic gas there are two basic concerns from 
a microscopic viewpoint. One, a thermodynamically consistent treatment of the 
coexistence of the two phases of nuclear matter (namely, the nuclear liquid and the 
surrounding gas) the other, a plausible description of the interface between 

the liquid and the gas []15| — pT7f1 . It has been shown that this problem can be taken 
care of by solving the density profile in the subtraction procedure of Bonche, Levit 
and Vautherin (BLV) [|T^,|T^]. In this method, the density profile of the liquid- 
plus-gas system and that of the gas are solved in a self-consistent procedure and 
then the extensive observables referring to the liquid are obtained as the difference 
between those of the two solutions. The influence of the surrounding gas on the 
surface energy is then automatically taken into account. The same concept can 
be applied to situations at zero temperature where drip nucleons occur and the 
system coexists with an outer nucleonic phase, as shown for very asymmetric cold 
semi-infinite nuclear matter in Ref. |L7| . 

In this paper we focus on the effect of the external gas on the structure and 
stability of finite nuclei with large neutron excess, to have a broader understand- 
ing of the conditions under which these nuclei may exist in neutron star matter. 
A preliminary study for isolated nuclei immersed in a nucleonic gas at zero tem- 
perature has been done recently [M; in the present work we extend these ideas in 
the appropriate astrophysical context. We also extend the calculations to non-zero 
temperature which is relevant at the formation stage of the neutron stars. 

We assume that in the nuclear matter at sub-nuclear density the nuclei are 
located in a lattice. In order to simplify our calculation, the Wigner-Seitz approx- 
imation || is applied, where each lattice volume is replaced by a spherical cell of 
radius R c , the nucleus being located at its centre. The matter in each cell is taken 
to be charge neutral, ie, the number of electrons is equal to the number of protons 
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in the cell. The neutrino density in this cell is determined from the /3-equilibrium 
condition. For a given average density of nuclear matter with a certain proton con- 
centration, the density distribution is solved self-consistently. The calculations are 
performed in the finite temperature Thomas-Fermi formalism. The matter in the 
cell in which the protons and neutrons coexist does not define the nucleus itself, 
the nucleus is identified after subtraction of the gas part generated self-consistently, 
as in previous investigations of excited nuclei or asymmetric semi-infinite matter 

In Section II, the model employed in the calculations is introduceed. The results 
and discussions are presented in Section III. Section IV contains the concluding 
remarks. 

II. MODEL 

In the inner region of the crust of a neutron star, we consider a mixture of 
neutrons, protons (some of which may exist as bound nuclear clusters), electrons 
and neutrinos in thermodynamic equilibrium. We contemplate both cold matter 
as well as matter at a finite temperature. We ignore the contributions from alpha 
particles and also from photons. We further ignore the plasma effects. All these 



effects are known to be rather small [211. The nuclear clusters are assumed to be 



arranged in a body-centered cubic lattice which we approximate by Wigner-Seitz 
(WS) cells defined as spheres with radius R c . Each cell is assumed to be electrically 
neutral and interactions among the cells are neglected. At the densities that we 
are interested in, the Fermi momenta of the electrons are much larger than their 
rest mass, the electrons are then extremely relativistic and can be assumed to be 
uniformly distributed in the cell. We assume the matter to be in /3-equilibrium, ie, 
the chemical potentials of the considered species fulfil 
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Pn = H p + He- (1) 

At a fixed temperature, under the condition of charge neutrality, the variables 
in the problem are the average baryon density (p), the proton fraction Y p — Y e (the 
electron fraction), the neutrino fraction Y v and the lattice radius R c . With the f3- 
equilibrium condition given by Eq. (|l|), if p n , p p and p e are known, then p v can be 
determined which in turn yields p v when the neutrinos are taken to be a degenerate 
Fermi gas. Then only three of the variables are left independent. The condition of 
/9-equilibrium of the baryons with the electrons and neutrinos guarantees that the 
system has the minimum free energy. If one imposes further the constraint that the 
lattice sites contain a particular nuclear species with a given charge Z c i and baryon 
number A c i (immersed in a nucleonic gas or not), then there is only one independent 
variable. 

To obtain the thermodynamic properties of the system of baryons, electrons 
and neutrinos, we minimize the free energy of the system in the WS cell with the 
constraint of number conservation of the individual species. Under the conditions 
of charge neutrality and /3-equilibrium, the relevant grand potential is given by 

n = F-j2v g A q , (2) 
<? 

where T is the free energy and the index q = (n,p) refers to neutrons and protons. 
The free energy has the following expression: 

F({p),Y p ,Y v ,T) = J [H(r) -Ts(r) + £ c (r) + f e (p e ) + f u (p„)]dr. (3) 

In Eq. (|3|) the integration is over the volume of the WS cell. Here 7i refers to the 
baryonic energy density excluding the Coulomb energy, s is the entropy density 
of the baryons, £ c the Coulomb energy density of the system and /„, f e are the 
free energy densities of the neutrinos and of the electrons (Coulomb interaction 
excluded) . 
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For the nuclear energy density, we use the Skyrme energy density functional. It 
is written as 
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where p = p n + P P and the effective mass of the nucleons m* is defined through 
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In the numerical calculations we shall employ the SKM* interaction, whose param- 
eters can be found in Ref. P^j . In the Thomas- Fermi approximation the kinetic 
energy density r q is given by 



r q (r) = -(S^pf 



at T = 0, 
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The fugacity r) q is obtained as 



J-3,/2(Vq) at T 7^ 0. 



(6) 
(7) 



77,(7-) = (p, - K(r))/T, 



(8) 



where V q is the single-particle potential experienced by nucleons (including the 
Coulomb part for the protons). The nucleonic density p q is related to r] q by 

3/2 

— I Ji/2(v q )- (9) 
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The functions Jk(Vg) i* 1 Eqs. (0) and @ are the standard Fermi integrals. 
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For the entropy density of the nucleons one has 

*(r) = EKV3) JsM/JiM - V,] p r (10) 
q 

The direct part of the Coulomb energy density £ c of the charged particles is 

£*(r) = \{ Pp {r) - Pe ) J ( Pp (r f ) - p e )^-^dv\ (11) 
while the exchange part is computed in the Slater approximation: 

= A *{P?{t)+pT)- (12) 
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For f e and f v we use the standard expressions [ 23| . For the minimisation pro 



cedure, we take recourse to the finite temperature Thomas-Fermi approximation 
which yields Eq. (|8|). This equation is solved self-consistently leading to the density 
distributions of neutrons and protons in the spherical WS cell. 

Under the conditions of interest here, the baryonic fluid in the WS cell can 
congregate into nuclei (we consider only the spherical shapes) located at the centre 
of the cell and may be embedded in a low density gas of nucleons. Near the edge of 
the cell, for hot or for isospin rich systems beyond the drip line, the baryon density 
profile is found to be practically uniform which is identified as the low density gas. 
To delineate the nucleus from the embedding low-density environment, we follow the 



BLV procedure which has been used earlier in the Hartee-Fock (HF) framework [18 



as well as in a Thomas- Fermi (TF) scheme for a hot nucleus in coexistence with 
the vapour surrounding it. The method is based on the existence of two solutions to 
the HF or TF equations, one corresponding to the liquid phase with the surrounding 
gas (LG) and the other corresponding to the gas (G) alone. 

For an isolated hot nucleus in equilibrium with the gas surrounding it, in the 
absence of Coulomb forces, the densities p q LG and the gas density p q G can be obtained 
from the variational equations 



it = (13) 

ins , , 

M = ° < 14 > 

where fi^ G and Q G are the grand potentials of the respective systems with the same 
chemical potential. In our procedure, the solution for p L = p LG — p G , which may be 
called the liquid profile, is independent of the box volume in which the calculation 
is done. In the presence of the Coulomb force, however, the Coulomb repulsion 
increases with the box size, driving the protons to the edge of the box and leading 
ultimately to a divergence problem. In order to work out a convergent prescription, 
Bonche et al. fl8| calculated the Coulomb energy E c from the subtracted proton 



density p P i G —p G aiJ -d obtained the density profile from the variation of the subtracted 
grand potential Cl = Q^g ~ ^g + E c with respect to both p q LG and p q G . 

In the astrophysical context, the situation is however different due to the presence 
of electrons. The Coulomb energies for the LG and G phases are 

E C LG = \ J (P P L G (r') ~ Pe)^^(p P LG (r) - Pe )dvdv\ (15) 
E G = \ J (PG(r) - Pe)^^(p P G (r) - Pe )dvdv' 

+ / Pli^j^jipUr) - Pe)dvdv'. (16) 

The direct part of the single-particle Coulomb potentials (8E1 G /8p p LG and 5E G /5p G ) 
for the LG are the G solutions is the same. It is given by 

V c d (r) = J (ff LG (r') - p^j-^dr'. (17) 

Since the system is charge neutral, the divergence problem does not arise. The 
solutions to the density profiles can be directly obtained from the variations of the 
total (with Coulomb) grand potentials Ql G for the liquid-plus-gas phase and Q G 
for the gas phase with respect to p\ G and p G , respectively. The resulting coupled 
equations are 
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Tvlcir) + Vl G + V£ G (p p LG , p e ) = N , (18) 
Tr f G ( r ) + V G + V<M G , Pg , p.) = av (19) 



At zero temperature they reduce to 

(37r2)2/3 2^ (PiG)2/3 + ^ + V ^LG,Pe) = ^, (20) 
(37r2)2/3 2^ ( ^ )2/3 + VS + V ^Plo,P P G,Pe) = AV (21) 

In Eqs. (|T8|)-(f2~l|), Vl G and V G refer to the nuclear part of the single-particle potential 
corresponding to the liquid-plus-gas and the gas solutions. V£ G (V G ) is the sum of 
the direct part of the Coulomb single-particle potential given by Eq. (|17]) and the 
exchange term — e 2 (3/7r) 1//3 (p^ G ^ G ,^) 1 / 3 for the protons in the LG (G) phase. At zero 
temperature, if the cluster is within the drip lines, the cell does not contain any 
nucleonic gas. Equations (|20|)-(|2lD then automatically lead to a gas solution which 
is zero throughout. 

The nucleonic chemical potential is determined from the constraint of nucleon 
number conservation. In the case of a given nuclear cluster with neutron number 
N c i and proton number Z c i embedded in a nucleonic gas, the conservation of the 
number of nucleons in the cluster requires [from Eqs. (|18D and ([19]) ] 

A** = J" {/ [Trilair) + V q LG {r) + V£q{t)) pl G {r)dv 

- J [T VG (r) + VS(r) + VS(r)} p 9 G (r)rfr} . (22) 

Here A q refers to N c i or Z c \. An equation similar to (|22|) follows from Eqs. (|20D and 
(pl|) at zero temperature. 

On the other hand, when (p) is given in a WS cell with a given proton concen- 
tration, the total number of neutrons iV and protons Z in the cell is defined. Then, 
the chemical potential is obtained as 
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^ = T I Mc(r) + vlG{r) + vlG{r)] ^(0*, (23) 



where now A q refers to iV or Z. 



III. RESULTS AND DISCUSSION 



In an earlier paper by some of the present authors |[20|| , calculations at zero 
temperature on some structural properties of isolated nuclei much beyond the drip 
line were reported. In that case, the excess pressure exerted by the enveloping gas 
stabilises the nucleus even beyond the nominal drip lines (which are defined by 
H n = or fi p = in the Thomas- Fermi approach). We now place these calculations 
in a broader context related to the environment existing in the inner crust of a 
neutron star. We take cognizance of the presence of electrons and neutrinos existing 
in /3-equilibrium with the neutrons and protons. We extend the calculations also 
at finite temperature. As mentioned earlier, the calculations are performed in the 
finite temperature Thomas- Fermi framework in a Wigner-Seitz lattice employing the 
SKM* interaction. For completeness, we also discuss the situation for the proton- 
rich nuclei though they may not exist in the crustal matter of the neutron stars. 

In the calculations on asymmetric infinite and semi-infinite nuclear matter in 



equilibrium with a drip phase, it was observed []l2" , |l6| , |r7| that the neutron-proton 
asymmetry could be increased arbitrarily till the two phases merge into a uniform 
system when the densities and the asymmetries of both phases become equal. On the 
contrary, in the self-consistent TF calculation for isolated finite nuclei, it was found 
that one cannot add or remove neutrons from nuclei arbitrarily PU| . There exists 
a limiting neutron-proton asymmetry / = (N c i — Z c {)/A c i beyond which the system 
becomes thermodynamically unstable; establishing chemical equilibrium between 
the nuclear phase and the gas phase further becomes impossible. 

In the realistic conditions considered in the present investigation, including the 
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lattice effects, we calculate the stability limits at T = as shown in the upper 
panel of Fig. 1. We find that both the limiting asymmetry (full line) and the 
drip asymmetry (dashed line) are not very sensitive to the atomic number Z c i of 
the nucleus. The influence of the lattice electrons on both the drip lines and the 
limiting lines for neutrons and protons is manifested in the lower panel of Fig. 1. 
The neutron lines are nearly unaffected. The proton lines are however influenced 
significantly, particularly for the heavier clusters. Both the proton drip and limiting 
lines get extended with inclusion of the lattice effects because of the dilution of the 
Coulomb force in the presence of electrons. The calculations are done in a lattice 
size R c = 15 fm. The results are found to be rather insensitive to the size of the 
lattice. 

The influence of temperature on the drip and limiting lines, taking into account 
the lattice effects, is displayed in Fig. 2. There, we compare the results obtained at 
T = and at T = 6 MeV. Both the neutron and proton drip lines are extended at 



finite temperature as found earlier in calculations for isolated nuclei p^ , |25| . On the 
other hand, the neutron limiting line shrinks with temperature, whereas the proton 
limiting line remains essentially unaffected. 

In Fig. 3 we present the density profile of neutrons (left panel) and protons 
(right panel) for an extremely neutron-rich nucleus 330 Pb (the neutron drip line is 
located at 276 Pb). We display calculations for both the case with lattice and the 
case without lattice (isolated nuclei) at T = 0. At a finite temperature T = 6 MeV 
we show the calculations only with lattice effects. The total density profile p\ G of 
the liquid-plus-gas phase is shown in the upper panels. For neutron-rich nuclei the 
influence of the electrons in the lattice is found to have only a nominal effect on the 
density distributions. The finite neutron density at the cell boundary even at T = 
reflects the presence of the neutron gas for a nucleus beyond the neutron drip line. 
With increasing temperature, the central density is depleted with the appearance 
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of a thickening tail. The gas densities p q G are shown in the bottom panels and the 
liquid densities p\ obtained as the difference between p\ G and p q G are shown in the 
central panels. 

The density profile p\ of the isolated cluster 330 Pb is found to be independent of 
the size of the box. However, in a lattice, because of the modification of the Coulomb 
force in the presence of the electrons, the nucleonic densities in the cluster are box- 
size dependent, though weakly. The liquid density vanishes at a distance r < R c 
even at T = 6 MeV. The neutron gas density is found to be practically constant 
throughout the box; at finite temperature, the gas density is larger as expected. 
For this neutron-rich system, at T = 0, the proton gas does not exist; however, at 
finite temperature (T = 6 MeV), a few protons are present in the gaseous state. In 
finite temperature calculations for isolated nuclei, the proton gas density profile is 
strongly polarised JT8 1 due to the repulsion from the nuclear core. In the presence 
of the lattice electrons, the density polarisation exists only in the vicinity of the 
core, at further distances the proton gas density is found to attain a nearly constant 
value. 

The density distributions for an extremely proton-rich isotope 140 Pb are shown 
in Fig. 4. The proton and neutron distributions at T = (both with and without 
lattice contributions) and at T = 6 MeV with lattice effects are displayed in the 
left and right panels, respectively. The electrons in the lattice dilute the Coulomb 
repulsion thereby lowering the proton chemical potential, particularly for proton- 
rich heavy nuclei, as an effect of which the proton drip line is extended. For instance, 
for an isolated Pb nucleus the proton drip line is located at 182 Pb; the influence of 
the lattice electrons pushes it to 143 Pb. The total (LG) proton density distribution 
for the isolated 140 Pb nucleus at zero temperature indicates the presence of a proton 
gaseous phase. This proton gas is strongly polarised as is clear from the proton 
gas density distribution shown in the bottom panel. With inclusion of the lattice 
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electrons, since the 140 Pb nucleus is just beyond the modified drip line, the proton 
gas density is extremely dilute as is evident from the bottom left panel. The central 
panel of Fig. 4 indicates that the liquid density profiles with and without lattice 
corrections at T = are not much different. The influence of the lattice electrons 
on the neutron density distribution is negligible as displayed in the right panel. For 
the present system, there is no neutron gas at zero temperature. The characteristics 
of the density distributions at finite temperature are found to be basically the same 
as discussed in the context of a neutron-rich nucleus in Fig. 3. 

All of the subsequent calculations we shall present are performed including the 
influence of the lattice electrons. In Fig. 5, the rms neutron and proton radii of the 
lead isotopes are displayed for T = as well as for T = 6 MeV. At zero temper- 
ature, except at the edges of limiting asymmetry, the rms radius changes almost 
linearly with the mass number. The faster change in radius near the edges points 
to the onset of instability. Similar observations were also noted in the calculations 
without lattice effects pO| |. Temperature increases the rms radii, its effect being 
more pronounced near the stability limit. By comparing the upper and lower panels 
of Fig. 5 (neutrons and protons) , one observes the growth of the neutron or proton 
skins with positive or negative asymmetry, the effect being more noticeable at zero 
than at finite temperature. 

We now explore the ambient conditions in which a particular nucleus can exist 
at a certain temperature in the inner crust of the neutron star. For this purpose, 
we have chosen 80 Ca and 170 Sn as representative systems. The nucleus 80 Ca is well 
beyond the Thomas-Fermi drip line ( 68 Ca being on the drip line) whereas 170 Sn is 
just at the drip boundary. As mentioned in Sec. II, among the four variables (p), Y p , 
Y v and R c , only one is independent once the N c i and Z c i values of the nucleus are 
constrained and /^-equilibrium assumed. For the existence of the particular isotope 
the values of the rest of the variables are then fixed. The correlations among (p), 
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Yp and R c are displayed in Fig. 6 for the aforementioned two isotopes at T = 0, 
3 and 6 MeV. The filled circles correspond to different values of R c at which the 
calculations have been performed. The values of R c range from 12 fm to 30 fm at a 
step of 3 fm, increasing with decreasing (p) for all the cases shown. The nucleus 80 Ca 
is beyond the neutron drip line, it is hence embedded in a neutron sea even at zero 
temperature. For the isolated nucleus, the density of the neutron gas is independent 
of the box size. Even in calculations with inclusion of effects from lattice electrons, 
the neutron chemical potential p n and hence Pq change very little with the cell size. 
This explains the gradual fall of (p) and Y p with increasing R c . For 170 Sn, since the 
nucleus is at the neutron drip boundary, the fall of (p) with increasing R c is vertical 
at T = as there is no neutron gas in the cell. At finite temperature, however, the 
cluster becomes embedded in a nucleonic gas of mostly neutrons and some protons 
causing a reduction in (p) and Y p with increasing R c . 

The ambient conditions for stability are investigated further for Pb isotopes 
ranging from A c i = 110 to A c \ = 330 at the same three temperatures mentioned 
for 80 Ca and 170 Sn. The calculations for Pb have been performed at a the value of 
cell size R c = 15 fm. The results are displayed in Fig. 7. The proton drip point for 
lead at T = is at A c i = 143, while the neutron drip point is at A c \ = 276. These 
points are indicated by the vertical arrows in the figure. In the upper panel of this 
figure, the average density (p) is plotted against mass number A d . The curve for the 
zero temperature case has a minimum at A d ~ 143 which is the proton drip point. 
Between the drip points (A d = 143 to 276), since there is no nucleonic gas around 
at T = 0, the average density (p) within the WS cell increases linearly with mass 
number. Beyond the drip points, the nuclei are surrounded by a nucleonic gas for 
their stability, (p) therefore increases faster. For a hot nucleus there are evaporated 
nucleons in the cell and the average density is larger. 

The proton fraction Y p as a function of the mass number is shown in the middle 
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panel of Fig. 7. A kink is seen to occur at A d = 143, the proton drip point below 
which Y p rises faster because of the presence of the proton gas. For A d > 143, the 
charge Z c i is fixed. Hence, with increase of mass number, Y p decreases linearly up 
to the neutron drip point beyond which the slope of the curve changes as neutrons 
appear in the drip phase. With temperature, because of evaporation of nucleons tak- 
ing place, the kinks smoothen out, but the overall qualitative features are the same. 
The lower values of Y p at finite temperature reflect the predominance of neutrons 
in the gas for isotopes that are not very proton-rich. The variation of the neutrino 
fraction Y v (= p v j (p)) with the isotope mass is displayed in the bottom panel. With 
increasing mass number, it is found to decrease. The effect of temperature on Y u 
is very weak. Inspection of Fig. 7 reveals that for a particular (p), A c i is generally 
double-valued. If the isotope is proton-rich, Y p is necessarily large and Y v also comes 
out to be large. For the same value of (p), the neutron-rich isotope is formed when 
Y v is small. 

So far our calculations have been done by fixing a given cluster (N c i, Z c i), which 
yields various sets of values for the parameters of the problem. A common practice 
is, however, to fix the parameters like (p) and Y p and study the change in the internal 
structure of neutron star matter. When /^-equilibrium is assumed, for a given value 
of (p) and Y p , the free energy density (equivalent to free energy per baryon in the 
present case) in the system of nucleons, electrons and neutrinos is a minimum in 
any given cell size R c . The cell size sets the periodicity in the neutron star matter; 
it is not a priori known and therefore is taken as a parameter. Even for given values 
of (p) and Y p , the neutrino concentration Y u (under /^-equilibrium conditions) would 
be different for different values of R c . One may look for a global minimum in the 
free energy density by varying Y v (through variation of R c ). This is displayed in 
Fig. 8 for T = and for three values of densities keeping the proton concentration 
Y p fixed at 0.3. The minima in the free energy density are marked by arrows. With 
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decreasing density, the cell size corresponding to the global minimum increases. 

To understand the origin of the minimum in the free energy profile presented in 
Fig. 8, the variation of its different components with lattice size for two values of 
(p) = 0.001 (left panel) and 0.01 fm -3 (right panel) are shown in Fig. 9. Here Y p is 
fixed at 0.3. Since these calculations are carried out at zero temperature, the total 
energy per baryon e tot = F/A is given by 

F/A = e N + e lat + e e + e v . (24) 

The different contributions ejy, e-iau e e and e u are the nuclear energy including 
the Coulomb interaction among the protons, the lattice energy (ie, the Coulomb 
interaction energy due to the presence of the electrons), the electron kinetic energy 
and the neutrino energy per baryon, respectively. It is seen that has a minimum 
at a lattice radius R c somewhat smaller compared to i?™ n , the value of R c at which 
F/A is minimum in Fig. 8 at the corresponding (p). This may be understood from 
the variation of the cluster size A c i with R c as shown in the middle panels of the 
figure. The cluster size increases monotonically with R c and for these neutron- 
rich clusters the minimum binding energy per nucleon occurs at around A c i ~ 100 
(for nuclei on the /3-stability line, this occurs for 56 Fe). The lattice energy e\ at 
decreases monotonically with R c . This is because for given Y p , the proton number 
and hence the electron number increases with increasing lattice size keeping the 
electron density constant, as a result of which e\ a t ~ —A 2 / 3 [BJ. The electron kinetic 
energy per baryon e e is a constant for a given density and Y p and therefore is not 
shown in the figure. The neutrino energy e u passes through a minimum at a value of 
R c somewhat smaller than R™ m where the neutrino fraction Y v is also a minimum 
(shown in the bottom panel). The competition between e\ at and the rising parts of 
(eN + e v ) determines the location of the minimum in the total free energy at T = 0. 
In Fig. 10, we display the thermal evolution of the cluster composition (Ad, Z c i) 
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at different fixed densities, at a particular value of Y p = 0.3 and for a lattice size 
R c = 15 fm. The number of neutrons and protons in the WS cell is then fixed. The 
full lines correspond to the mass number A c i of the cluster and the dotted lines refer 
to its charge. For the chosen conditions, the clusters formed at zero temperature 
in the WS cell comprise all the nucleons without any gas. As the temperature 
rises, the cluster size shrinks and becomes surrounded by the gas of the evaporated 
nucleons, the nuclear liquid and the gas being in thermodynamic equilibrium. The 
cluster evaporates completely at a particular temperature depending on the chosen 
density, and then the cell contains only the gas of nucleons and the leptons. We 
call it the boiling temperature T B . As the density (p) is decreased, the temperature 
Tb falls down. A similar situation was observed in the context of the liquid-gas 
phase transition in finite nuclei, where the phase transition temperature was found 
to decrease |26|] with an increase in the so-called freeze-out volume or a lowered 
average density (p). Figure 10 also tells us about the appearance of different nuclear 
clusters as the neutron star evolves in the formation stage. Initially, the temperature 
T may be as high as ~ 10 MeV and then all the nucleons will likely be in the gas 
phase. As time flows by, the system cools down and seeds of nuclei with increasing 
size start appearing from the gas. 



IV. CONCLUSIONS 



We have investigated the structural properties of nuclei and nuclear matter as 
can be found in the environment of the inner crust of a neutron star. These nuclei 
are dipped in a sea of electrons, neutrinos and of low-density nucleons. Because 
of the pressure exerted by the surrounding nucleonic gas, nuclei may exist even far 
beyond the nominal nuclear drip line. Following the BLV subtraction procedure, 
we give a prescription in a Thomas- Fermi framework to properly isolate the nucleus 
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from its environment. This is extremely helpful in further exploring the limits of 
stability. 

From our calculations, a limiting asymmetry in the neutron-protron concentra- 
tion emerges beyond which a nucleus even within the gaseous environment ceases to 
exist. A delicate balance between the Coulomb force and the diluted surface tension 
with increasing asymmetry and increasing density of the environment possibly plays 
a pivotal role here. We also investigate the ambient conditions, namely the average 
density, electron or proton fraction and the neutrino concentration, under which 
nuclei of a particular species can be formed at different temperatures in the stellar 
matter. We furthermore study the thermal evolution of the nuclear clusters at dif- 
ferent densities, which may serve as a guide to understand the nucleation process of 
different species of nuclear clusters from the nucleonic sea as the neutron star cools 
down in the earlier stages of its formation. 

We have left our calculations on a simple pedestal. That is, we have not included 
the extension of the Thomas-Fermi framework, we have not taken the shell-effects 
into account, possible plasma effects have been ignored and, similarly, the presence 
of a-particles at low densities has not been taken into consideration. We have worked 
with the SKM* interaction; at large asymmetry, its validity is not unquestionable. 
Sophistications in the approach or the use of more suitable interactions may change 
the results somewhat, but the qualitative features that emerge from our calculations, 
we believe, will remain mostly unchanged. 

Two of us (X.V. and M.C.) acknowledge financial support from the DGICYT 
(Spain) under grant PB98-1247 and from DGR (Catalonia) under grant 1998SGR- 
00011. 
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Figure Captions 

Fig. 1 The upper panel shows the drip line asymmetry and the limiting asymmetry 
as a function of the charge Z c i of the nuclear cluster in a lattice of size R c — 15 
fm at T = 0. The lower panel displays the same quantities in the N c i — Z c i 
plane with and without inclusion of lattice effects. 

Fig. 2 The drip and the limiting asymmetries in the N d — Z d plane at T = and 6 
MeV in a lattice of R c = 15 fm. The limiting proton asymmetry coincides at 
both temperatures. 

Fig. 3 The neutron (left panel) and proton (right panel) density profiles correspond- 
ing to the nuclear cluster 330 Pb are shown with inclusion of lattice effects at 
T = MeV (open circles joined by dotted lines) and at T = 6 MeV (dashed 
line). The full line corresponds to calculations for the isolated nucleus (with- 
out lattice) at T = 0. The upper panels refer to the liquid-plus-gas density 
(LG) and the central panels correspond to the liquid density (L). The lower 
panels show the gas density (G). 

Fig. 4 Same as in Fig. 3 for 140 Pb. Here the left panel corresponds to protons and 
the right panel refers to neutrons. 

Fig. 5 The neutron (top panel) and proton (bottom panel) rms radii at T = and 6 
MeV calculated in a lattice cell of size R c = 15 fm for Pb isotopes. 

Fig. 6 The average density and proton fraction of the nuclear matter in the Wigner- 
Seitz cell calculated for the nuclear clusters 80 Ca and 170 Sn. The size R c of the 
cell is varied from 12 fm to 30 fm at a step of 3 fm, as marked by the filled 
circles. The calculations are done for temperatures T = 0, 3 and 6 MeV. The 
full lines are drawn to guide the eye. 
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Fig. 7 The variation of the average density (p), the proton fraction Y p and the neu- 
trino fraction Y v as a function of the mass number A d of the Pb isotopes 
calculated in a Wigner-Seitz cell of size 15 fm. The calculations are done at 
T = 0, 3 and 6 MeV. The vertical arrows indicate the proton and neutron 
drip points. 

Fig. 8 The variation of the free energy per baryon F/A with the cell size at average 
densities (p) = 0.001, 0.005 and 0.01 fm -3 , the proton fraction being fixed at 
Y p = 0.3. The calculations correspond to T = 0. The vertical arrows refer to 
the minima in F/A. 

Fig. 9 The different components of the free energy per nucleon as a function of R c 
are shown in the upper panel for two different values of (p), at Y p — 0.3 and 
T = 0. The middle panels show the growth of the nuclear cluster with cell 
size. The bottom panels display the variation of Y v as a function of R c . The 
numbers on the left bottom panel are to be scaled down by a factor of two. 

Fig. 10 Thermal evolution of the cluster sizes for average densities (p) = 0.001, 0.005 
and 0.01 fm~ 3 . The full lines correspond to the mass number A c \ and the 
dotted lines refer to the atomic number Z c i. The calculations are performed 
with a fixed proton fraction Y p = 0.3 and a cell size R c = 15 fm. 
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